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ABSTRACT 

Using six-dimesional phase-space information from the Fourth Data release of the Radial Velocity Experiment (RAVE) over the range 
of Galactic longitude 240° < / < 360° and Vlsr < -239 kms -1 , we have computed orbits for 329 RAVE stars that were originally 
selected as chemically and kinematically related to oj Centauri. The orbits were integrated in a Milky-Way-like axisymmetric Galactic 
potential, ignoring the effects of the dynamical evolution of oj Centauri due to the tidal effects of the Galaxy disk on the cluster along 
time. We also ignored secular changes in the Milky Way potential over time. In a Monte Carlo scheme, and under the assumption that 
the stars may have been ejected with velocities greater than the escape velocity (V re i > V esc fi) from the cluster, we identified 15 stars 
as having close encounters with oj Centauri: (/) 8 stars with relative velocities V re i < 200 kms -1 may have been ejected ~ 200 Myr 
ago from oj Centauri; (ii) other group of 7 stars were identified with high relative velocity V re i > 200 kms -1 during close encounters, 
and seems unlikely that they have been ejected from oj Centauri. We also confirm the link between J131340.4-484714 as potential 
member of oj Centauri, and probably ejected ~ 2.0 Myr ago, with a relative velocity V re i ~ 80 kms -1 . 

Key words. Stars: abundances - Stars: kinematics and dynamics - Galaxy: structure - Globular clusters: individual (oj Centauri, NGC 
5139) 


1. Introduction 

Known as the most massive and luminous globular cluster of 
the Milky Way’s Halo, oj Centauri also shows unique physical 
properties in its structure, chemical enrichment and internal 
kinematics ( [Freeman & Rodgers ||1 975 [|Norris et al.||1997} |Lee| 


et al. 1999| Hughes & Wallerstein 2000{]Gnedin et al. 2002 


Bedinetal. 2004; So 


Mizut ani et al.||2003af Pancino et ai.||2003t |Hilker et al.l|2004 


2009 [ |Johnson & P ilachow ski|201 Of Bellini et al. 2010[ |Moriya 


lima et al. 2007, 2009a; Bellini et al. 


& Shigeyama||2010[ Da Costa |2012| ). The nature and origin 

of oj Centauri is still controversial. It has been proposed that 
oj Centauri is the remnant core of a tidally disrupted satellite 
galaxy, possibly destroyed long ago by the Milky Wa y (|Lee| 
et al. |1999[ |Majewski et al.] 2000; Be kki & Freeman||2003 ' 


Mizutani et al.||2003b |Ideta & Makino||2004[ |Tsuchiya et aL| 


2004[ |Romano et aL[2007 ; Sollim a et al.||2QQ9b| ), and some 

recent works suggest that it could have ori ginated from an even 
more massive system of about 10 9 M© (Valcarce & Catelan 

|20TT] ). 


those that have explored the surroundings of the cluster, looking 
for the controversial pair of tidal tails extending from it and 
reported by |Leon et aL| ( |2QQQ| , are the effects of differential 
reddening as spurious effects on the detection of tidal tails 
( [Law et aLl |2003 ), and the spectroscopy survey of 4 000 stars 
selected from the cluster red giant branch done by |Da Costa' 
& Coleman| ( |2008| ), and the more recent photometric survey of 
RR Lyrae stars over - 50 square degree s around the cluster by 
[Fernandez-Trincado et al.| ( [2013112015| ). None of these works 
found evidence of stellar overdensities beyond the tidal radius 
of oj Centauri. However, different results are inferred from the 
photometric analysis of the STREGA survey (Marconi et al. 
|2014| STRructure and Evolution of the GAlaxy), whose tracers 
are mainly variable stars (RR Lyraes and long-period variables) 
and main-sequence turn-off stars. These authors found evidence 
of stellar overdensities at - 1° from the cluster in the direction 
perpendicular to what was explored in the latter mentioned 
survey, nonetheless a more detailed chemical abundances and 
kinematics analysis would help to verify this result. 


Under this scenario it is expected that the orbit of such 
progenitor would have evolved along time because of the 
interaction with the Milky Way ([Tsuchiya et al.||2003| |Bekki 


& Freeman||2003| |Zhao||2004| |Ideta & Makino 2004| ) and part 
of its residual material would be distributed in the Galactic 
Halo. Several efforts have been made, with various degrees of 
success, to detect the relics of oj Centauri ’s progenitor. Among 


Stellar debris features associated with oj Centauri have also 
been identified recently in a wide range of Galactic latitude, 
and more specially in the solar neighborhood, with chemical 
patterns and ki nematics properties similar to inner populations 
of the cluster (Nissen & Schuster 2010} |Wylie-de Boer et al. 
|2010[ |Majewski et al.||2012| ), and more recently a group of 
RAVE stars with radial velocities and chemical abundances 
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consistent wit h stars observed fro m oj Centauri clus ter has been 
detected by ( [Kunder et al.||2014| ). As mentioned by |Meza et al.| 
(2005 ), the presence of nearby stars with odd kinematics and 
chemistry, suggest an extra-Galactic origin for them, possibly 
“tidal relicts”. 

From that starting point, we decided to investigate the kinemat- 
ical properties of stars having preferential retrograde motion 
(Vlsr < -239 kms -1 ), looking for the presence of co Centauri 


debris. We adopted an idea similar to that of [Pichardo et aL 


( |2012| ) and combined the accurate radial velocities, distances 
and [Fe/H] abundances from the Fourth Data Release of the 
Radial Velocity Experiment survey (RAVE DR4) by Kordo patis| 
et al.| ( |2013| ), with proper motions from the US Naval Obser¬ 


vatory CCD Astrograph Catalog (UCAC4) by |Zacharias et al.| 
( |2013| ), to select a sample of giant stars that are chemically and 
kinematically similar to oj Centauri. For each of star in the sam¬ 
ple, we integrate 10 5 orbits in a Milky-Way-like axisymmetric 
potential, to evaluate the probability of the occurrence of close 
previous encounters between the star and oj Centauri, within the 
cluster’s tidal radius. 



This paper is organized as follows: In §[2] we describe the 
data and methods for the selection of star candidates. In ^3] 
we briefly explore the effects of contamination in our sample 
by different Galactic components. In ^4] we describe the orbits 
model. Results and discussion are presented in §[5] Finally, 
conclusions are presented in ^6] 


2. Sample selection from RAVE 

The sample was selected from the RAVE DR4 catalog (Kor^ 
|dopatis et al. 112013] ), which provides accurate radial velocities 
with typical errors cr RV ~ 2 km s -1 , and distances and individual 
abundances with errors of about 10-20%, determined for 
approximately 390,000 relatively bright stars (9 mag < Idenis < 
13 mag). The proper motions on RAVE DR4 were compiled 
from several catalogs, however, in this work we use UCAC4 
( [Zacharias et al.|2013[ ). We used these data to make a kinemat- 
ical selection of RAVE stars possibly related to oj Centauri, 
taking also into account spatial distribution and metallicity, as 
well as some additional quality control cuts to select robust data. 


In this work we restricted our study to RAVE stars with 
Galactic longitudes 240° <1< 360°, where oj Centauri’s remnant 
candidates have been found (e.g., [Majewski et al.||2012| ). We 
restricted our sample to giant stars, with effective temperature 
between 4000-5500 K, and surface gravity 0.5 dex < log(g) < 
3.5 dex, following |Boeche et al.| ( [2011| ). Additionally, we 
required the stars to have high quality spectra (x 2 < 2000) with 
a signal-to-noise ratio S/N>20 ( algo_conv = 0 was required, 
indicating that the pipeline converges, see Kordop atis et al.| 
|2013| ). This cut allows us to obtain precise radial velocity 
measurements, typically cr RV < 2 kms -1 , in order to constraint 
the full space motion. The metallicity [Fe/H] distribution for 
giant stars within oj Centauri spans more than a magnitude order, 
from -2.2 dex < [Fe/H] < -0.7 dex ( [Johnson & Pilachowski| 
|2010| ), therefore we allow stars in our sample to be in this range 
of metallicity. 


A total of 9,024 RAVE stars in our final sample passed the 
selection and quality controls described above. Distances, ra¬ 
dial velocity and proper motions are all available for this sample. 


Fig. 1 . Toomre diagram for RAVE stars and their simulation with the 
Besangon Galaxy model. RAVE stars with quality control cuts (see 
are shown in grey dots, oj Centauri is shown with the green symbol. Iso¬ 
contours refers to Thin-Disk (green), Young Thick-Disk (black dashed), 
Old Thick-Disk (magenta) and Halo (cyan) simulated RAVE stars. Or¬ 
ange symbols (8 stars) with error bars correspond to RAVE stars having 
a high probability of encountering with oj Centauri and relative velocity 
V re i < 200 kms -1 (see ^5|, and red dots (314 stars) correspond to Halo 
stars having a negligible probability (< 0.01%) of encountering with oj 
Centauri (see & The vertical dashed line indicates zero Galactic rota¬ 
tion. (A color version of this figure is available in the online journal.) 


2.1. Kinematic selection 


With the RAVE kinematics and distances as input from the se¬ 
lected sample, galactic space velocities (U, V, W ) were computed 
according to the matrix equations of |Johnson & Soderblom| 
( |1987| ). For reference, we have adopted a right-handed co¬ 
ordinate system for (U,V,W), so that they are positive in the 
directions of the galactic center, galactic rotation, and north 
galactic pole, respectively. In this convention, we have adopted 
the solar motion with respect to the Local Standard of Rest 
( Uq,Vq,W q ) = (11.10,12.24,7.24) km s -1 , assumed that the 
Local Standard of Rest is on a circular orbit with a circular 
speed of 239 km s -1 , and the Sun is located at R Q = 8.3 kpc 
(e.g., |Brunthaler et al.|2011| ). 


Figure [I] shows the Toomre diagram for the Thin Disk, 
Young/Old Thick-Disk and Halo stars from the revised version 
of the Besangon Galaxy model (Robin et al.] |2014| ). The grey 
dots represent the 9,024 RAVE stars previously selected. As 
seen in the Toomre diagram, most RAVE stars are kinematically 
consistent with the different Galactic components. Then, stars 
with Vlsr < -239 kms -1 are more probably Galactic halo stars, 
and this work is focused on studying the orbit of these stars and 
their possible connection with oj Centauri. 


|Majewski et al.| ( |2012|) suggest that oj Centauri tidal debris 
are a primary contributor of retrograde velocity stars near the 
Sun. In fact, evidence has been found of field stars asso ciated 
with oj Centauri in the Galactic Disk by |Wylie-de B oer et al. 
(2010). The cut Vlsr < -239 kms -1 acts as kinematical mem¬ 
bership selection and helps us to obtain a sample of 329 stars 
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Fig. 2. RAVE stars in the Galactocentric radial velocity ( Vgrv ) vs. Galactic longitude space. Top panel refers to RAVE data, middle and bottom 
panel refers to Besangon Galaxy model-simulated RAVE data. a> Centauri is shown with the green cross. [Majewski et aL| ( |2012| ) stellar "streams" 
are plotted in blue symbols. Probable members found in this work are represented in orange symbols and are the same as in Figure [T] their 
size is proportional to the probability of close encounter with the cluster (see Sect. [4}, and RAVE Halo stars having a negligible probability of 
encountering with a> Centauri are shown with red dots. The black contours show the full RAVE sample and simulations, respectively. The coloured 
contours refers to the Galactic components from the simulation as in Figure [I] The grey scale refers to the logarithm of the number of stars in it, 
as indicated by the colour bar. (A color version of this figure is available in the online journal.) 


that is mostly free of contamination by the Galaxy’s Thin and 
Thick Disks. Given the magnitudes of these stars we are able 
to explore past dynamic links between these RAVE stars and co 
Centauri. For these stars we computed orbits using an axisym- 
metric model of the Galaxy scaling the model of Allen & Santil- 
Ian| ( |199l| (gravitational potential details can be seen at Q. 


3. Kinematic contaminants from the Milky Way 

In order to estimate the kinematical contamination in our sample 
by stars from the Galaxy, we simulated the RAVE data with the 
Besangon Galaxy Model. We used the new version described 
in [Robin et al.| p014[ ) where the Thick-Disk and Halo charac¬ 
terizations have been revised according to SDSS and 2MASS 
data. Moreover, the kinematics of the different populations is 
computed using a corrected asymmetric drift formula (Robin 
et al. 2015, in prep.) reliable out of the Galactic plane. The 
velocity ellipsoid of the Halo is assumed to be (131, 106, 85) 
km s -1 with the major axis pointing towards the Galactic center. 
To simulate the RAVE data, we run a simulation for all the 
directions and then a subsample of stars is randomly chosen in 
each magnitude bin in the I filter, containing exactly the same 
number of stars observed by RAVE. A full comparison of this 


new Besangon Galaxy Model with RAVE will be presented in a 
future paper (Robin et al. 2015, in prep.). 


In Figure [2j we present a comparison of the radial veloc¬ 
ity relative to the Galactic rest-frame (Vgrv) distribution for the 
simulated sample and the RAVE data. To obtain Vgrv , which 
is simply the star’s radial velocity as seen from the Sun’s loca¬ 
tion but in a reference frame at rest with the Galaxy, we use the 
formulation given by Hawkins et al. ( 2015| ) and the following 
equation: 


V G rv =RV+[U q cos (/) + (V© + V L y*) sin (/)] cos (b) + W 0 sin (b) 

( 1 ) 


The Besangon Galaxy Model reproduces very well the 
general kinematics observed in the RAVE data, as evidenced 
by the signal of the Disk observed in the sinusoidal-like density 
path seen in Figure |2j with a higher density around 1 = 210° and 
Vgrv —220 km/s (e.g., |Hawkins et al.|2015| ). 


Assuming a Gaussian velocity distributions and asymmetric 
drifts for the Thin-Disk, Young/Old Thick-Disk and Halo com¬ 
ponents, stars with V to tal >180 km s -1 have a high probability 
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of belonging to the Halo (see |Nissen & Sc huster 2010, and 
references therein). However, in the case of the Besangon 
Galaxy model ( [Robin et al.||2014| ), the velocity distribution of 
the Thick-Disk is non-Gaussian with an extended tail toward 
Vtotal < 239 kms -1 , then in our simulated RAVE stars, those 
with 180 kms -1 < Vtotal < 239 kms -1 might belong to the 
Thick-Disk population. In this sense, since our selected stars 
have Vlsr < -239 km s -1 , they are most probably free of 
Thick-Disk contamination (and certainly free of the Thin-Disk 
one), nonetheless we could expect Halo stars with preferential 
retrograde motion in the Galactic longitude 240° < l < 360° 
explored in this work. 


4. Orbits in an axisymmetric Galactic model 

In this work, we ran 10 5 pairs of orbits for each RAVE star 
and co Centauri clusterjpair, using the Monte Carlo schem^] 
describe in Pichar do~et al.| ( |2QT2] ). For each Monte Carlo set of 
orbits pairs, the proper motions, radial velocities and distances, 
with their respective error bars, are used for both the star and 
the cluster. As initial test, we calculate orbits backward in 
time over 1 Gyr from the current position and velocity of co 
Centauri in a Milky-Way-like axisymmetric Galactic potential 
(dynamical friction, spiral arms and bar effects are ignored), 
which is the | Allen & San tillan (1997]) model (see details in 
|Pichardo et al.| T012| ), scaled with the new values Rn = 8.3 


Brunthaler et al. 


(2011). Each 


kpc, 0o = 239 kms -1 given by 
Monte Carlo run employs a different scaled Galactic potential, 
using the uncertainties in Rq and 0o. Then, we calculated the 
probability of close encounters in the past, which was defined 
as: (0 having a minimum approach distance, d m i n < 100 pc (e.g., 
[Ferndndez-Trincado et al.|20T5] ), such that during the encounter 
the RAVE star is within the region of bounded stellar motions 
of co Centauri, (ii) having an encounter time, r enc <0.2 Gyr; in 
which the bar shape of the Galactic bulge and its effect on the 
Galactic orbits may be approximately ignored, this assumption 
is supported by the results in Appendix [A] (We verify that the 
introduction of a bar in our potential, using the prolate bar 
model of |Pichardo et al.| |2004| , does not change significantly 
our conclusions for eight of our candidates with relative ve¬ 
locities V re i < 200 kms -1 listed in Table [2J) . to Centauri has an 
orbital period r or b ~ 0. 08 - 0.12 Gyr (e.g7 |Dinescu et al.|1999 


Miz utani et al.||200~3b] ) thus, in 0.2 Gyr the cluster will have 
approximately two perigalactic points. To illustrate the orbital 
behavior in about 0.2 Gyr, Figure [3] shows an example of a pair 
of orbits star-cluster for a RAVE star having a close encounter 
(d m in <100 pc) with co Centauri in r enc ~ 0.2 Gyr and a relative 
velocity of V re i = 53.6 kms -1 . The average distance star-cluster 
is, d avg = 0.9 kpc, from t- 0 Gyr until the close encounter time 
(f ~ Tenc\ 


For co Centauri itself, we calculate the orbit, based on the 
current distance of 5.57 ± 0.08 kpc (Navarrete et a l. 2015| ), proper 
motion ji a cos(S ) = -5.08 + 0.35 mas yr -1 and us = -3.57 + 0.34 


mas yr 1 , radial velocity 232.2+0.7 km s 


-l 


1999), and position a =13h26m47.24s and 6 = 


(Dinescu et al. 
47°28 , 46.5" 


1 Uncertainties are considered as lcr variations, and a Gaussian 
Monte Carlo sampling generates the parameters to compute the 
present-day positions and velocities of the co Centauri and RAVE stars 

2 The probability is defined as: Prob = Number of orbits having close 
encounters with the cluster / N to taU where N tota i - 1 x 10 5 Monte Carlo 
simulations 
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Fig. 3. Orbit projection in the potential axisymmetric for 
"J1 31340.4-484714" (blue line) and co Centauri (red line) in X, 
Y Galactic coordinates on the Galactic plane. The solar position is 
(X Q , Y 0 , Z 0 ) = (8.3, 0, 0) kpc in t = 0 Gyr. The solid curves show the 
orbits integrated backwards from initial position in t - 0 Gyr (open 
symbols) to r enc = 0.2 Gyr (filled symbols) having a close encounter 
with d m i n <100 pc. (A color version of this figure is available in the 
online journal.) 


( |Harris|1996| ). 


We briefly highlight some limitations of our calculation, 
since we ignore the effects of: ( i ) dynamical evolution of the 
cluster, and tidal effects on the cluster exerted by the Galactic 
Disk in each orbit; (ii) mass loss, as co Centauri was likely 
more massive in the past (see [Bekki & Freeman|2003[ [jdeta & 
Makino 2004 ), this reduces the probabilities we find. We have 
also ignored secular changes in the Milky Way potential over 
time ( [Lind et al.|2015| ). 


5. Results and discussion 

Within the sample of 329 stars with Vlsr < -239 kms -1 , a 
group of 15 RAVE stars were identified as having a probable en¬ 
counter with co Centauri within 0.2 Gyr at 100 pc distance, and 
another group of 314 stars as having unprobable (< 0.01%) close 
encounter with co Centauri with the same distance and time limit. 

For those 15 stars with probable encounters, we identified 7 
stars with high relative velocities V re i > 200 kms -1 , for which 
it is unlikely that they have been ejected from co Centauri; 
and 8 stars were identified with relative velocities V re i < 200 
kms -1 , which are more likely associated with co Centauri 
(both groups of stars are listed Table [2]). These stars (8 more 
likely) have an average Galactic rotation velocity of Vlsr ~ 
-280.9 + 37.5 kms -1 , which is close to that of the co Centauri, 
Vlsr = -265.3 + 55.0 kms -1 , as computed by us from the 
cluster’s data. The resulting values Ulsr,Vlsr , and Wlsr 
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for these stars are given in Table [5] Table [5] gives the same 
parameters for the 314 stars classified as unlikely to have been 
ejected from oj Centauri. 


We identified 314 RAVE stars with unlikely (or null) 
encounters with oj Centauri occurred within 100 pc. Figure [2] 
shows that their Galactocentric radial velocity (red symbols) is 
in agreement with the kinematics predicted by the Besan 9 on 
Galaxy model for the Halo. Their [Fe/H] and Vlsr are in 
agreement with characteristics of the inner Galactic halo (e.g., 
[Dinescu et al.|1999) . 


Figu re |2| s hows stellar debris candidates (blue symbols) 
from [Majewski et al.| ( [2012] ), with a kinematical and chemical 
pattern linked to oj Centauri. In the same figure, stars in our 
kinematical and chemical RAVE selection (orange symbols) 
have retrograde orbits similar to that of oj Centauri, and are 
kinematically coherent with Majewski et al7|( |2012 ) sample. It is 
also evident that our 15 RAVE stars candidates with retrograde 
velocities do not follow the Disk kinematics in both Figures [T] 
and [2 


We also identified in our sample one star with 
RAVEID=J131340.4-484714, previously identified beyond 
the cluster tidal radius by |Anguiano et al.| ( |2015| ) and listed as 
potential member of the oj Centauri globular cluster. Under 
our numerical simulation scheme, we found 3,582/100,000 
orbits for J131340.4-484714 having close encounters with 
oj Centauri, i.e. our results suggest that the star has been 


ejected from oj Centauri with an ejection velocity (V e 


ejection 


80 


kms -1 ). Figure [5] shows that this star could have been ejected 
approximately ~2.0 Myr ago, following an orbit very close to oj 
Centauri as is evident in the Figure [6] . This peculiar star could 
also be a rare "no-escaper" i.e. a star on a stable orbit outside 
the cluster tidal radius ( [Ross et al.|l997| . 


5.1. Close encounters: Star-Cluster 


Figures [5} and [6] shows the relative velocity (V re l) distribution 
for 15 RAVE stars having close encounters in the past, and the 
average distance ( d avg ) distribution for the pairs of orbits from 
t - 0 Gyr to t - r enc , respectively. As is evident in Figure [5j 
the vast majority of the close encounters are produced during 
approximately time intervals less than two times the orbital 
period of oj Centauri, with relative velocities V re i < 200 kms -1 
(e.g., [Lind et al7| |2015| ) in most cases, and closed orbits with 
small values of d avg , such that the pairs of orbits (star - cluster) 
travel with close distances between them (see ilustrated example 
in Figure [3]), i.e., this approach clearly suggests that their orbits 
are confined to the oj Centauri orbit. Table [ 2 ] lists 15 RAVE stars 
with the more frequent relative velocity ( V re i ), average distance 
(d aV g), and enconter time (r enc ) for each Monte Carlo simulation. 


The escape velocity of oj Centauri cluster at the center is 
about V esc ,o = 60.4 km s -1 , and at the cluster half-mass radius 
V esc ,h = 44 km s -1 (e.g., 


Gnedin et al. 2002). Given the obser¬ 


vational errors, and relative velocity V re i > V esc ,o star-cluster 
encounters, it is seems unlikely that the star has been tidally 
stripped from the cluster, as this mechanism is not efficient to 
explain such fast encounters in our simulations. For V re i > V esc $, 
the star most possibly has been ejected from the cluster, which 
is explained in some scenarios by binary system interactions 
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Notes. 1 Heliocentric radial velocity ( RV ). Tablets published in its entirety in the electronic edition of the journal. A portion is shown here for guidance regarding its content. 
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Fig. 4. 2D density plot for [X/H] and Vgrv » in the Gala ctic longitude range 240° < / < 360°. The orange and red symbols are the same as in Figure 
[T] The cyan dots refers oj Centauri members from |Johnson & Pilachowski] ( [2010) . (A color version of this figure is available in the online journal.) 


and black holes (e.g. |Pichardo et al.|2012[ |Lind et al.|2015 
references therein). 


see 


5.2. Chemical distribution for probable oj Centauri members 

Figure [4] shows Vgrv vs. [Al/Fe], [Si/Fe], [Ti/Fe], [Ni/Fe], 
[Fe/H] and [Mg/Fe], for the 15 RAVE stars (orange symbols) 
and the whole RAVE sample (in grey bins), The distribution for 
the 15 stars found as most likely debris from the oj Centauri pro¬ 
genitor, clearly follows the oj Centauri chemical patterns, as ob¬ 
served for member_stars within oj Centauri ’s tidal radius (e.g., 
Johnson & Pilachowski|2010[ and references therein). The heavy 
elements (Si and Ti) generally follow the typical enhancement of 
~ +0.3 dex found in the cluster population, oj Centauri exhibits a 
complex history of chemical enrichment, with an abundant metal 
poor and intermediate population (e.g, 

|2010| ). The extended enhancement in c 
served in our sample of 329 RAVE stars is not conclusive to 
select oj Centauri members from RAVE data, however, we have 
shown on the Toomre diagram of Figure[l] that these stars follow 
the Halo kinematics. It should be noted that oj Centauri spans 
ranges in all the analysed abundance patterns which are similar 
to those of the Milky Way, so any similarity between these quan¬ 
tities is not (poorly) useful to confirm the origin of these stars. 


Johnson & Pilachowski 
lemical abundances ob- 


5.3. Statistical significance 

To establish the statistical significance of our results we compare 
the Monte Carlo simulations (see 0 to five mock samples 
created using the Besangon Galaxy model (Robin et al.| |2014| ). 
For each mock sample we chose 329 stars randomly satisfying 
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Table 3. Mock sample from Besangon Galaxy model (Robin et al. 
|20T4] >. 


Id 

Number of false positives (Nfalse) 

Mock sample 1 

4 stars 

Mock sample 2 

2 stars 

Mock sample 3 

2 stars 

Mock sample 4 

4 stars 

Mock sample 5 

1 stars 


the selection criteria defined in §(2j and the RAVE/UCAC4 
selection bias (errors in radial velocities, proper motions and 
distances are added). The number of false positive are listed in 
Tabl43] 


We consider a high significance in our detection if 


U~significance ~ Ntrue 


{Nfalse) > 3 X CFf a i se 


( 2 ) 


where, N true = 15 (see text), and as we have 5 mock sam¬ 
ples from Table [3j we take the average over all the samples 
((Nfalse) = 2.6), obtaining a mean and standard deviation 
(o'false = 1.2) for this values. Then, we find a significance 
crsignificance = 10.33. Therefore, the number of stars reported 
in the present paper are cr significance > 3xcr fa i se compared to 
the mock sample. We can conlude that the methodology in §[4] 
robustly show that there is a statistically significant population 
of stars from RAVE survey possibly ejected from oj Centauri. 
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Table 2. Derived ejection parameters for the 15 RAVE stars more likely having close encounters with oj Centauri. 


RAVEID Axisymmetric Potential Non-axisymmetric Potential 



Vrei 

d a vg 

Tenc 

Vrei 

d a vg 

' r enc 


kms -1 

kpc 

Myr 

kms -1 

kpc 

Myr 

Encounters with V re i < 200 kms 1 

J131340.4-484714 

75.9 

0.25 

2.4 

79.7 

0.26 

2.1 

J120407.6-073227 

170.3 

2.25 

37.4 

143.2 

2.28 

36.1 

J151658.3—123519 

137.9 

2.3 

35.4 

139.5 

1.88 

32.1 

J143500.9-264337 

108.3 

1.32 

71.4 

108.2 

1.35 

71.1 

J143024.9-085046 

105.6 

2.34 

28.4 

182.7 

2.34 

28.1 

J144618.8-711740 

166.8 

2.89 

88.4 

132.9 

2.87 

97.1 

J144734.5-722018 

77.0 

1.34 

45.4 

187.6 

1.81 

43.1 

J123113.7-194441 

167.8 

2.31 

30.4 

170.4 

2.31 

31.1 


Encounters with V re i > 200 kms 1 


J150703.5-112235 

224.7 

2.36 

23.4 

227.9 

2.39 

24.1 

J040133.8-832428 

225.2 

2.34 

29.4 

237.7 

2.34 

29.1 

J144926.4-211745 

200.8 

2.34 

39.4 

229.9 

2.79 

39.0 

J153016.8-214028 

678.0 

2.84 

118.4 

718.6 

3.48 

122.1 

J152554.5-032741 

206.0 

2.84 

48.4 

323.7 

4.02 

57.1 

J124722.9-282260 

207.3 

2.38 

46.4 

208.2 

2.85 

48.1 

J110842.1-715260 

275.0 

1.34 

45.4 

291.5 

1.91 

45.1 


6. Conclusion 

We identified 8 RAVE stars at galactic longitudes 
(240° < / < 360°), with [Fe/H] metallicities in the range 
-2.2 dex <[Fe/H]< -0.7 dex, with Galactic retrograde motion 
similar to that of oj Centauri, which could have had a close 
encounter with oj Centauri with a minimum approach d m i n <100 
pc during the last 200 Myr in the past. Our numerical scheme 
suggest that these stars have been probably ejected from the 
cluster with a velocity V ejection < 200 kms -1 , and their kinemat¬ 
ics connection with oj Centauri makes them good candidates to 
be part of the predicted sequence of tidal streams in retrograde 
rotation left behind the cluster over time. 

The probabilities we estimated may neglect some physical 
evolution of the orbit and/or internal structure of the oj Cen¬ 
tauri progenitor, but the agreement with other results in the 
literature supports our conclusion that the 8 RAVE stars re¬ 
ported in this study are candidates to be oj Centauri stellar debris. 

We have also identified another group of 7 stars having close 
encounters with oj Centauri, however the high relative velocities 
V re i > 200 kms -1 during the encounters, makes it unlikely that 
they have been ejected from oj Centauri. 

The other 314 stars with retrograde orbits have a negligible 
probability of close encounters within 100 pc in 0.2 Gyr. 
We have classified this group of stars as less likely to have 
ties with oj Centauri. Their chemical distribution and Galactic 
retrograde rotation are similar to those of the inner Galactic halo. 

We have not found an evident external structure beyond 
the tidal radius of the cluster from RAVE data on Galactic 
retrograde rotation. 

In order to facilitate the reproducibility and reuse of our 
results, we have made available all the orbit simulations in a 
public repository 


3 http: //£ ernandez-trincado.github.io/ 
Fernandez-Trincado/simulations.html 
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Fig. 5. Distribution of relative velocities (V re i) during close encounters for the whole integration time, t = 1 Gyr (pink histograms), and t < 0.2 Gyr 
(black histograms). (A color version of this figure is available in the online journal.) 
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Fig. 6. d aV g distributions, defined as the average distance between pairs of orbits through the Galaxy from the present-time (or t - 0 Gyr) until the 
close encounter time (r enc ). The pink distribution refers to the whole integration time (1 Gyr), the black distribution refer to close encounters with 
relative velocities within 0.2 Gyr and 100 pc distance. (A color version of this figure is available in the online journal) 
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Appendix A: Orbits in a non-axisymmetric Galactic 
model 

We also run the simulations in a non-axisymmetric galactic po¬ 
tential, in order to verify the bar effect in our results. The Galac¬ 
tic model that we considered is the analytical potential presented 
in §[4| and the non-axisymmetric Galactic potential including 
a prolate bar, given by Pichardo et al.| ( |2004| to the galactic 
bar, with a total mass of ~1.6xl O i0 and angular v elocity 
of Q b = 55 ±5 kms _1 kpc _1 (e.g., [Pichardo et al. 2012). Using 
the same considerations that ^4] to model the orbits of the RAVE 
stars and co Centauri, 10 5 pairs of orbits were integrated to 0.2 
Gyr in the past. Figure |A 1| and Figure |A.2| shows the results to 
the stars listed in Table |2| the probabilities of close encounters 
are similar in both Galactic potential, and the bar effects for in¬ 
tegration time <0.2 Gyr, does not affect the results presented in 
this work. 
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Fig. A.l. Relative velocity (V re i) distribution for 10 5 pairs of orbits integrated backward in time over 0.2 Gyr in an axisymmetric (pink histograms) 
and non-axisymmetric potential (black histograms). 



d avg [kpc] d avg [kpc] d avg [kpc] d avg [kpc] d avg [kpc] 


Fig. A.2. d avg distribution. The pink and red histograms are the same as in Figure [Ad] 
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